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We discuss the relation between remnant fragmentation in inelastic high-energy hadronic inter- 
actions and muon production in extensive cosmic ray air showers. Using a newly developed tool, a 
simple and flexible hadronic event generator, we analyze the forward region of hadronic interactions. 
We show that measurements of the Feynman-a; distribution in the beam fragmentation region at 
LHCf will be key to understanding muon production in air showers quantitatively. 
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I. INTRODUCTION 



Muon production in cosmic-ray induced air showers 
has been of great interest recently. The number of muons 
observed by the Auger collaboration exceeds the one 
predicted by popular hadronic interaction models by as 
much as 50% 

A possible mechanism for more muon production in 
air-shower models has been proposed by the EPOS 
model which describes successfully a large variety 
of RHIC data. In that model, the production of most 
hadron species (in particular, of baryons) relative to neu- 
tral pions is enhanced. Thus, a larger fraction of the 
shower energy remains in the hadronic channel which, in 
turn, leads to a higher number of muons than in other 
hadronic air shower models such as QGSjet-II and 
SibyllQ. 

In this paper we point to another mechanism that in- 
fluences muon production in air showers significantly, 
namely the treatment of the remnant of the projectile 
in a hadronic interaction. The break-up of the hadronic 
remnant determines the amount of energy remaining in 
the hadronic channel of the air shower (by producing 
baryons or charged mesons) relative to that "lost" to the 
electromagnetic part due to production of neutral pions. 
Remnant break-up shows up in the very forward region 
(large Feynman-xj?) of collider experiments and therefore 
measurements of the LHCf collaboration Q will be im- 
portant to understand muon production in air showers. 

A generic difficulty in studying hadronic interactions 
for air showers is the fact that one always compares inher- 
ently different Monte-Carlo event generators. As typical 
result of detailed comparisons one one has model A which 
produces e.g. more muons than model B. One can then 
compare generic properties of the underlying hadronic 
models, like multiplicity, baryon production, strangeness 
or charm production, forward scattering and try to argue 
for an observed effect. This method remains always very 
qualitative, since one cannot disentangle the importance 
of the different contributions and one is always left with 
the uncertainty that other (possibly neglected) effects are 
actually also important. 

In order to relate a specific property of air showers (e.g. 
muon number) unambiguously to a specific feature of a 



given hadronic interaction model (e.g. remnant break-up) 
we developed a new hadronic interaction model. The goal 
is to keep the model as simple and flexible as possible. 
We then study air shower properties by changing only one 
parameter at a time. To investigate the influence of the 
remnant break-up on muon production we explore differ- 
ent treatments of the remnant but do not vary anything 
else. Likewise, the effects from baryon production can be 
studied by exclusively varying the diquark probability in 
the string fragmentation. 



II. THE TOOL, THE PQCD EVENT 
GENERATOR PICCO 

The new tool, PiCCO (pQCD Interaction Code for 
COsmics), is a standard event generator. Details will 
be published elsewhere, here we just outline the basic 
features. The interaction is modeled by a superposition 
of a soft and a semi-hard Pomeron. The Pomeron pa- 
rameters are fit to reproduce the total and elastic cross 
section for hadron-hadron scattering. Hard scattering is 
performed by Pythia and the partons are mapped 
onto strings. Each Pomeron gives exactly 2 strings each 
of which connects a quark-antiquark pair scattered out 
of the interacting particles. The longitudinal momentum 
distribution of the string-ends is the same as in other 
models, dxj^fx for each (anti)-quark. The remaining 
particle, the remnant, has a x^'^dx distribution in the 
case of baryons, versus x~^'^dx for mesons Q. Fragmen- 
tation of the strings is done via the Lund fragmentation 
scheme as implemented in Pythia. 

The treatment of the remnant is important for air 
shower applications. With a given probability Pox, the 
remnant can be in an excited state with an invariant 
mass M; the distribution for M is dhP/M"^. To model 
the fragmentation of such an excited remnant state, it 
could be viewed, for example, as a quark matter droplet 
or as a simple string. In our model we choose the string 
approach. If the remnant has a remaining light-cone mo- 
mentum fraction x'^ (which is the momentum fraction of 
one string-end) , the corresponding momentum fraction of 
the other string-end is given by x~ = M^/a;"*". Remnant 
strings may have different fragmentation parameters in 
order to enhance forward baryon or strangeness produc- 
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FIG. 1: (Color online) Feynman-x distributions of protons 
for proton-proton collisions. The data is from Ref. The 
model featuring extreme (100%) remnant excitation probabil- 
ity, and the diquark approach of Pythia, lead to a dip in the 
forward spectrum. 



tion. 

The remnant treatment outlined above is not 
new, it has been used in many models: QGSjet- 
98/01 [|, QGSjet-II g, the versions NeXus 2 [| and 
neXus 3 and EPOS Q. In other models, the rem- 
nant is modelled by attaching one of the strings result- 
ing from Pomeron exchange to a diquark. While this 
approach has some drawbacks which will be outlined be- 
low, it is also theoretically less appealing, because one of 
the interactions (the one including the diquark) is pre- 
ferred over the other ones. By separating the remnant, 
all interactions are treated on equal footing. 

III. SOFT REMNANT BREAK-UP 

In this section we consider soft remnant break-up, i.e. 
remnant excitation. We shall first discuss the relevant 
parameters. Pex is the excitation probability. In the case 
of excitation, the invariant mass is determined randomly 
according to the distribution dJVP /AP inbctwccn the lim- 



FIG. 2: (Color online) Feynman-x distributions of tv" for tt^- 
proton collisions. The data is from Ref. [l^ 



its A/iiiin and Mmax- In our approach, wc choose M^i^ to 
be the minimum mass corresponding to the flavor content 
of the remnant plus twice the pion mass. In the case of 
a proton- like remnant (quark content uud), for example, 
Mmin = (0.938 + 0.28) GeV. M^ax is taken to be 0.5^/s. 
Our simulations show that the mass limits are more of a 
technical nature, their exact value does not change the 
results significantly. 

The excitation probability pox, however, is of utmost 
importance. The exact value of this parameter has been 
fitted to reproduce low-energy forward scattering data 
from Ref. see Fig.[TJ We find Pex = 0.6. Also shown 
in the same figure is the forward particle distribution for 
the case Pox = 1, i-e., when the remnant is excited always. 
Notice how this affects the forward distribution of final- 
state protons: a dip appears in the region 0.5 < xp < 1. 
The explanation is simple. If the remnant decays into 
two or more particles, the region below xp < 0.5 will 
be populated. The region 0.5 < xp < ^ corresponds 
to a mere longitudinal momentum loss of the remnant 
(since it does not contain the "participant" quarks from 
the incoming hadron) rather than to its fragmentation. 
The lower panel of Fig. [1] depicts the forward particle 
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distribution at higher energies. 

Pythia implements a different approach. Instead of 
treating the remnant separately, one of the strings cre- 
ated in the interaction is attached to a diquark from the 
projectile while the other string is attached to the re- 
maining valence quark. This topology, however, appears 
to have problems explaining the ratio of anti-Omega to 
Omega baryons. Sufficiently massive strings ending in a 
diquark tend to produce more anti-baryons than baryons 
near that end, which is not supported by the data p^ . 
In Fig. [1] we see that a diquark string also fails to repro- 
duce the forward proton production data, since, similar 
to a decaying remnant, the momentum of the diquark is 
shared among two or more particles. 

For pion induced reactions, the situation is different, 
with the remnant consisting of a simple quark-antiquark 
pair. Again, with a probability pcx, the remnant gets ex- 
cited and decays by forming a string. A tt'^ +p reaction at 
250 GcV has been studied by the EHS/NA22 collabora- 
tion [l^. Fig. [2]shows the forward production of charged 
particles. In the case of PiCCO, a better fit is obtained 
with a lower excitation probability. Epos, which has a 
higher excitation probability yields a even slightly bet- 
ter description of the data. The forward tt^ production 
is shown in the lower panel of Fig. [2l Here, enhanced 
breakup improves the description of the differential cross 
section, but all models are still far from the data. In fact, 
no other model is able to describe this data. If one mod- 
ifies parameters extremely to reproduce neutral pions, 
then one fails to describe the charged pion spectra. The 
situation is similar with data on tt^ proton scattering at 
360 GeV measured by the NA27 LEBC/EHS Collabora- 
tion [l^ . A new measurement of production in pion 
induced reactions would be extremely helpful. 



IV. INFLUENCE ON MUON PRODUCTION IN 
AIR SHOWERS 

How does the remnant affect muon production in air 
showers? The basic reasoning is very simple: if the rem- 
nant decays into two or more particles its energy is re- 
distributed among particles which stay in the hadronic 
channel and eventually produce muons (charged mesons, 
baryons) and particles which go into the electromagnetic 
channel (mainly neutral pions). 

Let us first consider the break-up of baryons. The rem- 
nant of a baryon induced collision is always a baryon. 
The flavor content might of course change due to charge 
exchange (for example, a proton induced reaction might 
give a neutron as remnant, which nevertheless remains 
in the hadronic channel of the air shower). If the rem- 
nant however gets excited and breaks up, generically it 
produces a number of neutral pions. Therefore, baryonic 
remnant break-up should in fact reduce muon produc- 
tion. 

For charged-pion induced reactions, the outcome is 
very different. Due to charge exchange, the remnant 



can be either of the flavor of a charged pion (ignoring 
strangeness for the purpose of this argument) or of a neu- 
tral pion. The decay of an excited charged pion therefore 
reduces muon production since additional neutral pions 
are produced. The opposite is true for the decay of neu- 
tral pion remnants, where additional charged pions in- 
crease the energy available in the hadronic channel. The 
important question is which effect is stronger: do charged 
pion remnants produce more tt^ than neutral pion rem- 
nants produce tt* ? We start by considering low multi- 
plicities, when the decaying remnant produces only two 
particles. With one Pomeron exchange, the remnant is 
with 50% probability charged and with 50% probability 
neutral. This can be easily seen by inserting a quark- 
antiquark pairs (underlined) into a ud system: 

projectile remnant -I- stringcnds 

ud ^ ud + dd ( 1 ) 

— *■ uu + ud (2) 

— > du + uu (3) 

dd + du (4) 

(HI) and ^ will give a charged remnant, ^ and ^ a 
neutral one. The decaying remnant will try to produce 
an equal amount of charged (tt"*" , 7r~ ) and neutral mesons 
(tt*^, 77 ''■'), since again a break-up will produce with the 
same probability a or a dd-pair. But the production 
of T]*-'' is suppressed due to its large mass, therefore we 
will find more charged particles in the final state than 
neutral mesons. The consideration of strangeness and 
baryon production enhances this effect. For high multi- 
plicities we know that we get 1 /3 neutral and 2/3 charged 
pions, therefore the fraction of tt° is strongly reduced. 
We see that in total, the overall number of neutral pions 
is reduced due to break-up. Therefore, this mechanism 
will, in all, enhance muon production. Any additional 
production of strange mesons or of baryons will amplify 
this effect in favor of particles which stay in the hadronic 
channel of an air shower. 

Fig. [3] shows the energy fraction of neutral pions in the 
remnant, for the two extreme cases without excitation 
(Pex = 0) and with complete break-up {pex = !)■ We 
see that remnant break-up indeed reduces the energy in 
neutral pions. 

Next, we calculate the total muon number with a full 
air shower simulation using Seneca[i^. We simulate 
proton induced average vertical showers, the observation 
altitude is 1400 m above sea level, and we normalize to 
the QGSjet-II model. The result is shown in Fig. 51 
for the different cases of enhanced break-up {pex = 1) 
of baryonic and mcsonic remnants. We observe that ad- 
ditional break-up of baryons reduces muon production 
while break-up of mesons enhances it. The net effect is 
an enhancement of muons, since pion induced collisions 
are more abundant. 

We are now in a position to consider as well enhanced 
baryon production. The parameters for the diquark 
production in the string fragmentation are adjusted to 
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FIG. 3: (Color online) Energy fraction of neutral pions in the 
remnant for ■K'^+Air collisions, for the cases with and without 
break-up. Remnant break-up reduces the energy fraction in 
Tv" and thus enhances muon production in air showers. 
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FIG. 5: (Color online) The ratio of total muon numbers of ver- 
tical proton induced air showers for different remnant break- 
up and baryon production scenarios. Enhanced baryon pro- 
duction means pdiq.str = 0.12 and pdiq.rom = 0.3 instead of the 
default value pdiq = 0.10 from Pythia. 
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FIG. 4: (Color online) The ratio of total muon numbers of 
vertical proton induced air showers for enhanced break-up 
(Pex = 1.0 instead of pcx = 0.6) of baryonic and mesonic 
remnants. 



give qualitatively similar results to EPOS. We choose 
Pdiq,str = 0.12 and J3diq,rcm = 0.3 as diquark-anti-diquark 
pair production probabilities in central strings and rem- 
nant strings, respectively. The higher probability for di- 
quark production from remnant strings is motivated by 
the remnant encountering a denser target and leads to 
more baryons in the forward region, as found to be im- 
portant by the authors of EPOS. Fig. [5] shows that 
muon production in PiCCO with default parameters (i.e. 
Pdiq.str/rom ~ 0-1 ffom Pythia) is within 10% similar 
to QGSjet-II. Enhanced baryon production gives 25% 
more muons. However, the combined effects of both en- 
hanced baryon production and remnant break-up add up 
to 40%. 

At low energies, complete remnant break-up of baryons 



is excluded by the data, (see Fig. [T]) and the results 
shown in Fig. [5] should actually be considered as an up- 
per limit. However, we do not know whether the rel- 
evant parameter changes at high energies. If we keep 
the one obtained at low energies, Pox = 0.6, we find a 
rather flat distribution xpdn/ dxp of protons at 2 TeV 
center of mass energy, whereas one notices a dip in the 
forward scattering spectrum for the case of complete rem- 
nant break-up. 



A motivation for an enhanced remnant break-up is the 
fact that at higher energies the projectile probes smaller 
gluon momenta in the target and therefore encounters a 
higher gluon density. This effect has already been in- 
vestigated within the Color Glass Condensate (CGC) 
framework for hadron-nucleus collisions at colliders in 
Ref. [l^ and has been applied to air showers in Ref . [13 . 
The main consequence of the enhanced remnant break- 
up is a suppression of forward particle production; this 
leads to a faster absorption in air showers and hence to 
a lower shower maximum Xmax- Furthermore, an en- 
hanced muon production was observed in Ref. but 
was attributed mainly to an increased overall multiplic- 
ity. Efforts are currently ongoing to implement this par- 
ticular mechanism of forward suppression into the PiCCO 
model, which will allow us to test for its influence on the 
muon number. We will also investigate projectile break- 
up due to a dense target in proton-proton collisions at 
LHC energies, and apply this mechanism to air shower 
simulations. 
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This experiment plans to measure photons, neutral pions 
and neutrons in the very forward region with a small de- 
tector between the two beam pipes of the LHC collider. 
Except for the diffractive peak, neutrons and protons arc 
expected to show similar distributions in this region, with 
neutrons being somewhat outnumbered. The results are 
shown in Fig. [S] The slope of the forward neutron spec- 
trum will give insight into remnant break-up at these 
energies and will therefore allow one to constrain muon 
production in air showers. It is also interesting to exam- 
ine the absolute normalization of the very forward neu- 
tron spectrum as this relates directly to the treatment of 
the flavor content of the remnant. QGSjet-II predicts 
less neutrons since it allows for only one charge exchange 
reaction. 

As for neutral pions: here the differences due to rem- 
nant break-up are much smaller. The relatively flat for- 
ward spectrum of EPOS stems from a specific imple- 
mentation of the popcorn effect [l^ , where the two lead- 
ing particles of string fragmentation can be interchanged 
with some probability. 
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We have shown that remnant fragmentation in 
hadronic interaction models for air showers has signif- 
icant influence on muon production. If the projectile 
remnant decays into two or more particles the energy 
fraction which goes into the hadronic or electromagnetic 
channel changes. Enhanced baryon break-up reduces and 
enhanced meson break-up increases muon production in 
air showers. At low energies, forward scattering data 
constrains the break-up probability of the remnant to 
roughly 60%. The measurements of the LHCf experi- 
ment will be extremely helpful for our understanding of 
remnant break-up at higher energies, and hence will pro- 
vide important constraints for muon production in air 
showers. 



FIG. 6: (Color online) Feynman-x distributions of identified 
particles for proton-proton collisions at 14 TeV cms energy. 



V. FORWARD SCATTERING AT LHC 
ENERGIES 

Finally, we compare model predictions in the kinematic 
range relevant for the forthcoming LHCf Q experiment. 
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